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Abstract 

We report all-electron, density-functional calculations with large Gaussian polariza- 
tion basis set of the recently synthesized octahedral O B24N24 cage that is perfectly 
round by symmetry, and boron-nitride (BN) clusters that its existence might sug- 
gest. We consider whether it is energetically possible that the two halves of this 
round cage could cap the BN nanotubes, modeled by B28N28 and B32N32 . The 
energetics show that BN nanotubes with such round caps, are only slightly less 
favorable than the BN clusters containing six squares as the only defects in the 
otherwise perfect hexagonal lattice. A larger BggNgg O cage formed from B24N24 by 
adding sufficient hexagons to isolate all squares is not very favorable energetically. 
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Boron-nitride clusters have been made by vacuum arc melting, dissolved in 
pyridine, and detected using time-of-flight mass spectroscopy, which shows 
that the B24N24 cluster is made in abundance [1]. Its structure is likely to 
be the largest, and thus most round, alternate boron-nitride (BN) cage in 
which all atom pairs are equivalent by octahedral symmetry, in analogy to the 
icosahedral Cgo cluster. The latter cluster is a fullerene, which is a hexagonal 
graphite sheet made to close into a cage by twelve pentagons as required by 
Euler's theorem [2], for which there are no alternate BN cages [3]. Alternate 
BN cages are made by substituting alternate sp'^ carbon atoms of certain 
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fuUerene-like cages with boron and nitrogen atoms. These cages contain only 
squares, hexagons, octagons, etc., and only strong BN nearest-neighbor bonds. 
The simplest and most stable alternate BN cages contain hexagons that are 
closed by exactly six isolated squares [4]. However, B24N24 is the alternate 
cage for the perfectly round C48 cluster with Oh symmetry [5] . Both are made 
round by the addition of twelve defects-six octagons and the corresponding 
six squares required by Euler's theorem. When made from carbon this cluster 
is quite unfavorable [6], but apparently not when made from boron nitride. 
Our calculations show that the B24N24 cage containing six squares isolated by 
hexagons is about 4-5 eV lower than the octahedral B24N24 cage[4]. Perhaps 
BN cages are almost as versatile as fuUerene cages; that could be the case if 
defects in boron nitride cages cost significantly less energy than they do in 
fuUerene cages. 

A study of the smaller fuUerenes suggests that it is never favorable to replace, 
in the context of Euler's theorem, two pentagons with a single square [7]. 
So in carbon cages squares are likely to be twice as energetically unfavorable 
as pentagons, which cost, when isolated, at least 2.3 eV each [8]. Heptagons 
and octagon are less energetically favorable in carbon cages, but are required 
to bend and twist free-standing nanotubes [8] . The energetics of fuUerenes if 
applied to B24N24 would argue strongly against its existence. 

The most abundant fuUerenes are buckminsterfuUerene, Ceo, and C70. The 
latter can be viewed as being formed by cutting buckminsterfuUerene into 
hemispheres in a highest-symmetry direction and inserting a ring of ten atoms, 
and then rotating the top hemisphere half a revolution before reconnecting it. 
It is natural to think about extending that process to larger capsules [10] and 
ultimately to the (5,5) carbon nanotube in the infinite limit [12,11]. Contrary 
to early belief, Cqo is not the most stable fuUerene. The larger icosahedral 
fuUerenes are more stable [13], but Ceo is the most stable in a wide range of 
cluster sizes and shapes [14] . Perhaps perfect roundness facilitates making Ceo 
and B24N24 . 

If, on the other hand, extra defects in boron nitride fuUerenes are so costly 
that they never occur except in B24N24, then boron nitride caps have a very 
special shape. If half a defect-free BN fuUerene caps a BN nanotube just as 
half a defect-free fuUerene caps a carbon nanotube [15] then each BN cap 
has precisely three squares. Furthermore, the three squares must lie on the 
periphery of the nanotube, and, as three points determine a plane, the BN 
nanotube cap must be flat. Defect-free flat and triangular BN caps are quite 
different from hemispherical carbon nanotube caps. 

If B24N24 is spherical and under some circumstances abundant, it is natural to 
ask if an additional ring of BN atoms can be added to form a less abundant but 
still possible cluster that is reminiscent of C70, and if further rings are possible 
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to make hemispherically capped BN nanotubes. We study the energetics of 
this possibihty using analytic density functional theory. 

Our calculations use the Slater-Roothaan method. It uses Gaussian bases to 
fit the orbitals and Kohn-Sham potential [16] of density-functional theory. 
The method through robust and variational fitting is analytic and variational 
in all basis sets [17]. The most general functionals that it can treat so far 
are certain variants of the Xa functional [18]. In particular, it can handle 
different ct's on different elements analytically and variationally so that the 
atomized energies of any cluster can be recovered exactly, and all energies 
are accurate through first order in any changes to any basis set. Thus the 
energy does not change if we choose to constrain any fit, which we do for 
all fits at no cost to the calculated energetics and negligible cost in com- 
puter time. Our work has spawned research efforts that brought both use- 
ful fitting techniques and density-functional theory to Gaussian-based ab ini- 
tio quantum chemistry. Thus there are good basis sets for fitting the Kohn- 
Sham potential. The best s-type fitting bases are still those scaled from the 
s-part of the orbital basis [19]. A package of basis sets has been optimized 
[20] for use with DGauss [21]. We use the valence double-C orbital basis set 
DZVP2 and the pd part of the (4,3;4,3) (A2) charge density fitting basis. 
Ahlrichs' group has generated a RI-J basis for fitting the charge density of 
a valence triple-C orbital basis set used in the Turbomole program [22]. We 
use its pdf part with the standard 631 IG** basis set downloaded from PNNL 
(http://www.emsl.pnl.gov:2080/forms/basisform.html). For both sets of cal- 
culations the values of a chosen for B and N are, respectively, 0.786751 and 

0. 767473, which give the exact atomic energies [23] when used with the triple- 
C/RI-J basis. The accuracy of this method using Hartree-Fock values of a [24] 
away from the sides of the periodic table is roughly that of Becke's half-and- 
half method [25], but as it requires no numerical integration. The method is 
ideal for use with Gaussian basis sets [17]. Forces are computed nonrecursively 
using the 4-j generalized Gaunt coefficients [26]. 

The symmetry restricted optimization of the BN cages was performed using 
the BFGS computer code [27]. The optimized hemispherically capped BN 
nanotubes are shown in Fig. 1 while their symmetries are given in the Table 

1. The (4,4) nanotubes are generated by extending along the four-fold axis 
[Fig 1(b)]. The caps start off perfectly round by the symmetry of B24N24 . 
They remain quite round after geometry optimization. In the figure they look 
a little fiat at the tops and bottoms because the octagons take up a significant 
portion of their caps. The (3,3) BN nanotube generated by extending the 
three-fold direction would have rounder caps, but are probably less likely to be 
made because they contain significantly more curvature. The binding energy 
per BN pair (BE), the energy gap between the highest occupied molecular 
orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), and 
the vertical ionization potential calculated at the optimized geometries are 
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Fig. 1. (Color online) The optimized structures for hemispherically-capped BN 
nanotubes. (a) and (b) are two different views of the B24N24 cage: (a) along the C3 
axis, (b) along C4 axis, (c) B28N28 (C^h) cage obtained by adding a ring of eight 
alternating B and N atoms, (d) tubular B32N32 (Ss) cage by inserting two rings of 
eight alternating B and N atoms (see text for more details). 

also given in Table I. 

The coordinates (in A) of the B and N atoms optimized with 6311G**/RI-J 
basis set for the recently synthesized octahedral B24N24 cage are (0.66, 1.66, 
2.59) and (2.74, 1.65, 0.74) respectively for the three-fold axis in the (1,1,1) 
direction. The positions of all remaining of atoms in the cage can be generated 
by symmetry operations. The HOMO and LUMO both have e symmetry and 
yield a large gap of 4.94 eV that agrees with the value of 4.95 eV [1] obtained 
using DV-Xa at the PM5 geometry. 

The C4/1 B28N28 cage [Fig. 1(c)] is generated from the base B24N24 cage by 
cutting the latter into two halves after orienting it along the C4 axis, then 
inserting a ring of eight alternating B and N atoms perpendicular to the axis, 



Table 1 

The binding energy per BN pair (BE), the energy gap between the highest occupied 
molecular orbital and the lowest unoccupied molecular orbital, and the vertical 
ionization potential (VIP) for the optimized BN cages, all in eV. 



Symmetry 6311G**/RI-J DZVP2/A2 
BE GAP VIP BE 



B24N24 





14.63 


4.94 


6.66 


14.64 


B28N28 


Cih 


14.89 


4.86 


8.48 


14.90 


B28N28 


T 


14.95 


5.38 


8.64 


14.95 


B32N32 


Ss 


14.96 


4.95 


8.36 


14.98 


B36N36 


Td 


15.06 


5.33 


8.42 


15.11 


BgeNge 





14.97 


4.25 


7.74 


15.12 




(^) (b) 

Fig. 2. The optimized structures for the favorable B28N28 (left) and B36N36 (right) 
cages. 

i.e. horizontally, and then rotating the top half by an eighth of a revolution. 
The resultant cage structure has C4/1 symmetry with eight inequivalent atoms. 
If another ring of four BN dimers is inserted in the same manner, then the 
resultant B32N32 cage is a tubular structure with Sg symmetry [Fig. 1(d)]. 
Table I shows systematic increase in the binding energy by going from B24N24 
to the B28N28 cage (0.26 eV per BN pair) and from B28N28 to B32N32 cage 
(0.06 eV/BN pair). The successive additions of more BN rings will result in 
further stabilization of the BN tubular cages and will ultimately lead to the 
(4,4) BN nanotube. Thus our total energy calculations do not rule out the 
possibility of BN tubes with round caps similar to those observed in case of 
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Fig. 3. (Color online) The optimized structures for the (a) B24N24 and (b) BggNgg 
cages. 

carbon nanotubes. A recent molecular dynamics study of the growth mecha- 
nism of BN nanotubes has also noted the formation of a cap containing four 
squares and octagon for this nanotube [28]. Molecular models of caps with 
four squares, octagons and hexagons were also speculated earlier [30]. All the 
BN cages studied have rather large HOMO-LUMO gaps of about 4.9-5.5 eV 
similar to those observed for the infinite BN nanotubes (~ 5.5 eV)[29]. The 
vertical ionization potential (VIP) calculated as the energy difference of the 
two self-consistent calculations for the neutral and singly positively charged 
cage also show relatively large values. 

Fig. 2 shows two examples of the most favorable BN cages that are composed 
of eight equivalent triangular faces meeting at six isolated squares [4,3]. Their 
BE, HOMO-LUMO gap and VIP are given in Table 1. The smaller cage, T 
B28N28 , Fig. 2(a) has the same stoichiometry as the middle-sized nanotube 
of Fig. 1. That C4/1 B28N28 hemispherically-capped nanotube is related to 
this cage by the insertion of twelve defects: six squares and six octagons. 
Using the triple-^ basis the energetic cost is 0.14 eV per defect and using the 
double-C basis the cost is 0.12 eV per defect, which is more than an order of 
magnitude less than the cost of any additional defect in a carbon fullerene. 
The calculations indicate that these cages are energetically more stable than 
the B24N24 cage. The trend is similar to that observed in case of Ceo where 
higher fuUerenes have larger binding energy than the C6o[13]. 
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Fig. 4. (Color online) The possible round cap of the (8,8) nanotube based on the 
BgeNge cage. 

Cgo does not cap the most important carbon nanotubes. Probably the metal- 
hc (10,10) nanotube is most important [31]. It can be capped by halves of 
round, but not perfectly round like Ceo, icosahedral C240 [32]. B24N24 can be 
enlarged by adding hexagons. Another such round cage of octahedral sym- 
metry is BgeNgg . The two optimized round cages of this study are shown in 
Fig. 3. The Table shows that BggNgg is the more stable than the B24N24 . Its 
halves (Cf. Fig. 4) would cap a round (8,8) BN nanotube. It is clearly differ- 
ent, from B24N24 , in that while being mostly round, its twelve squares stick 
out significantly, like the detonators of a sea mine. Furthermore, BggNgg is 
not substantially more stable than the other cages of Table I, despite the fact 
that it is at least three times as large. Thus we think it unlikely that hemi- 
spherical caps for larger than BN (4,4) nanotubes, if they exist at all, will 
have isolated squares, but no really attractive possibilities come to mind. For 
example, six square-octagon-square units could be placed on the two-fold axes 
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of the octahedral group, but the Euler overclosure would have to be balanced 
by two extra octagons inserted elsewhere to form a hemisphere. Of course, BN 
nanotubes could be capped by hemispheres of carbon fuUerenes. 

To summarize, all electron analytic density functional based calculations have 
been performed on selected boron-nitride cage structures including the re- 
cently synthesized B2/1N2/1 cage. The energetics of the BN tubular structures 
with the end caps based on this cage points to possible existence of round 
capped BN nanotubes. Amongst the BN cages that satisfy the six isolated 
square rule the higher B28N28 and B32N32 cages are found to be energetically 
more stable than the smallest B24N24 consistent with earlier calculations [4], 
but apparently they are not more abundant experimentally [1]. Additional 
topological defects in the most favorable BN cages seem more energetically fa- 
vorable than defects in the fuUerenes. Probably round caps on BN nanotubes 
with circular cross section are more likely than fiat caps on triangular carbon 
nanotubes. 

The Office of Naval Research, directly and through the Naval Research Labo- 
ratory, and the Department of Defense's High Performance Computing Mod- 
ernization Program, through the Common High Performance Computing Soft- 
ware Support Initiative Project MBD-5, supported this work. 

References 

[1] T. Oku, A. Nishiwaki, I. Narita, and M. Gonda, Chem. Phys. Lett. 380 (2003) 
620. 

[2] H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl, and R. E. Smalley, Nature 
318 (1985) 162. 

[3] H.-Y. Zhu, T. G. Schmalz, and D. J. Klein, Int. J. Quantum Chem. 63 (1997) 
393. 

[4] M.-L. Sun, Z. Slanina, and S.-L. Lcc, Chem. Phys. Lett. 233 (1995) 279. 

[5] A. D. J. Haymet, Chem. Phys. Lett. 122 (1985) 421. 

[6] B. I. Dunlap and R. Taylor, J. Phys. Chem. 98 (1994) 11018. 

[7] P. W. Fowler, T. Heine, D. E. Manolopoulos, D. Mitchell, G. Orlandi, R. 
Schmidt, G. Seifert, and F. Zerbetto, J. Phys. Chem. 100 (1996) 6984. 

[8] B. I. Dunlap, Int. J. Quantum Chem. 64 (1997) 193. 

[9] B. I. Dunlap, Phys. Rev. B 50 (1994) 8134. 

[10] R. C. Haddon, L. E. Brus and K. Raghavachari, Chem. Phys. Lett. 131 (1986) 
165. 



8 



[11] S. lijima, Nature, 354 (1991) 56. 

[12] J. W. Mintmire, B. 1. Dunlap, and C. T. White, Phys. Rev. Lett. 68 (1992) 
631. 

[13] B. I. Dunlap, D. W. Brenner, J. W. Mintmire, R. C. Mowrey, and C. T. White, 
J. Phys. Chem. 95 (1991) 8737. 

[14] T. G. Schmalz, W. A. Seitz, D. J. Klein, and G. E. Kite, J. Am. Chem. Soc. 
110 (1988) 1113. 

[15] M. Fujita, R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 45 
(1992) 13834. 

[16] W. Kohn and L. J. Sham, Phys. Rev. 140 (1965) A1133. 
[17] B. I. Dunlap, J. Phys. Chem. 107 (2003) 10082. 

[18] E. C. Vauthier, A. Cosse-Barbi, M. Blain, and S. Fliszar, J. Mol. Struct. 
(Theochem) 492 (1999) 113. 

[19] B. I. Dunlap, J. W. D. Connolly, and J. R. Sabin, J. Chem. Phys. 71 (1979) 
3396; 4993. 

[20] N. Godbout, D. R. Salahub, J. Andzelm, E. Wimmer, Can. J. Chem. 70 (1992) 
560. 

[21] J. Andzelm, E. Wimmer, Physica B, 172, 307 (1991); J. Chem. Phys. 96 (1992) 
1280. 

[22] K. Eichkorn, F. Weigend, O. Treutler, and R. Ahlrichs, Theor. Chem. Acc. 97 
(1997) 119. 

[23] A. Veillard and E. Clementi, J. Chem. Phys. 49 (1968) 2415. 
[24] K. Schwarz, Phys. Rev. B 5 (1972) 2466. 
[25] A. D. Bcckc, J. Chem. Phys. 98 (1993) 1372. 
[26] B. I. Dunlap, Phys. Rev. A 66 (2002) 032502. 

[27] W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, 
Numerical Recipes The Art of Scientific Computing (Cambridge University 
Press: Cambridge, England, 1986) p 309. 

[28] X. Blase, A. De Vita, J.-C. Charlier, and R. Car, Phys. Rev. Lett. 80 (1998) 
1666. 

[29] X. Blase, A. Rubio, S. G. Louie, and M. L. Cohen, Europhys. Lett. 28 (1994) 
355. 

[30] M. Terrones, W. J. Hsu, H. Terrones, J. P. Zhang, S. Ramos, J. P. Hare, R. 
Castillo, K. Prassides, A. K. Cheetham, H. W. Kroto, and D. R. M. Walton, 
Chem. Phys. Lett. 259 (1996) 568. 



9 



[31] A. Thcss, R. Lcc, P. Nikolaev, H. Dai, P. Petit, J. Robert, C. Xu, Y. H. Lee, 
S. G. Kin, A. G. Rinzler, D. T. Colbert, G. E. Scuseria, D. Tomanek, J. E. 
Fischer, and R. E. Smalley, Science 273 (1996) 483. 

[32] J. A. Harrison, S. J. Stuart, D. H. Robertson, and C. T. White, J. Phys. Chem. 
B, 101 (1997) 9682. 



10 



